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Abstract. Pulsational properties of 1.8 M© stellar models covering the latest stages of contraction toward the 
main sequence up to early hydrogen burning phases are investigated by means of linear nonadiabatic analyses. 
Results confirm that pre-main sequence stars (pms) which cross the classical instability strip on their way toward 
the main sequence are pulsationally unstable with respect to the classical opacity mechanisms. For both pms and 
main sequence types of models in the lower part of the instability strip, the unstable frequency range is found 
to be roughly the same. Some non-radial unstable modes are very sensitive to the deep internal structure of the 
star. It is shown that discrimination between pms and main sequence stages is possible using differences in their 
oscillation frequency distributions in the low frequency range. 
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1. Introduction 

In the vicinity of the main sequence, the HR diagram is 
confusing as stars of similar global properties but different 
stages of evolution lie at the same position. In some cases 
the young pre-main sequence stars are recognized through 
specific characteristics, for instance the presence of nebu- 
losity or high degree of activity. But, this is not always the 
case. An alternative is to take advantage of seismological 
pieces of information whenever it is possible. We therefore 
carry out a comparative study of the pulsational behavior 
of pre-main sequence (hereafter pms) and main sequence 
(hereafter postZams) stars of intermediate mass, in the 
classical instability strip. 

Our study is restricted to the latest stages of the pre- 
main sequence contraction when traces of the formation 
phases have already disappeared, so that a pms model can 
be built with the quasi-static approximation. 

The outer layers of pms and postZams stars having al- 
most the same effective temperature and gravity are very 
similar. As these layers drive the pulsation in this range of 
temperature, it is reasonable to expect that pms stars in 
the instability strip are also destabilized by classical opac- 
ity mechanisms and that the same type of modes than for 
the postZams stars are excited. This idea has been con- 



firmed by Suran (1995 ) and in a more detailed work but for 
pure radial modes of very young and cool stars by Marconi 
& Palla (199S). This supports the vibrational instability 
as an explanation for the observed periodic variations de- 
tected for a few stars which are suspected to be in the 
pre- main s equence stage (Kurtz & Marang 1994 , Catala 
et al. 



1997, Kurtz, 1999). 
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Differences, on the other hand, exist between the two 
stages: whereas pms has some relicts of its gravitational 
contraction phase, the main sequence nuclear burning 
modifies the inner core structure, and postZams mod- 
els develop c hemica l inhomogeneities. In a preliminary 
work, Suran ( |l998|) has stressed some consequences of 
these differences in the structures of the frequency distri- 
butions. The present paper further investigates and com- 
pares the pulsational properties of these two phases for 
non-radial oscillations. It is organised as follows: stellar 
modelling and its physical inputs are described in Sect. 
2, where comparable pms and postZams models are se- 
lected. Oscillation frequencies associated with nonradial 
modes for both types of objects are compared in Sect. 3 
and 4, and emphasis is put on the structure of the fre- 
quency spectra in both stages. Finally in Sect. 5, we dis- 
cuss the possibility to infer the evolutionary stage, i.e. the 
age of a variable object in this region of the HR diagram, 
through the pattern of its frequency spectrum. 
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2. Static models and evolutionary sequences 

2.1. Evolutionary tracks 

In the early stages of pms evolution, various rapid- dynam- 
ical and thermal- processes occur which modify the inter- 
nal structure of the protostar. By the time the protostar 
reaches the classical instability strip on its way toward the 
main sequence, these complex phenomena have long disap- 
peared and the evolution has considerably slowed down. 
The usual quasi-static approximation then holds and is 
assumed here. Rotation is not considered and the star is 
modelled in a spherically symmetric equilibrium state. 

Equ ilibri um models have been computed with CESAM 
(Morel 1997 ), an evolutionary code with a numerical ac- 
curacy of first order in time and third order in space. The 
models are built with the fast version of the code i.e. 300 
mesh points in mass; pulsation calculations on the other 
hand use an extended 2000 mesh grid, calculated from the 
static model by means of the basis of the spline solution. 



The EFF equation of state (Eggleton et al. [L973J) 
is used as it is sufficient for our purpose in this mass 
range. Opacities are provided by the Livermore Library 



(Iglesias & Rogers 1996), complemented at low tempera- 
tures (T < 10000K) by Alexander & Ferguson 's ( |l994|) 
tables. The nuclear network contains the following species 
1 H, 3 He, 4 He, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0. The main nu- 
clear reactions of pp+CNO cycles are included with the 
species 2 H, 7 Li, 7 Be set at equilibrium. Nuc lear r eaction 
rates are taken from Caughlan & Fowler ( 198§| ); weak 



screening is assumed. The isotopic helium ratio is fixed 
at He 3 /He 4 = 1.01 10 -4 and t he mi xture of heavy ele- 
ments is solar (Grevesse & Noels 1993[ ). In the convection 
zones, the temperature gradient is computed according to 
the standard mixing-length theory, with the mixing-length 
defined as I = aH p , where H p is the pressure scale height. 
Assuming a constant value of a along the main sequence 
(Fernandes et al. 1998), we have chosen a = 1.67 , the 
value obtained for a calibrated Sun computed with the 
same physical assumptions. The atmospheric layers are 
treated with the Eddington's grey approximation. 

We consider the lower part of the instability strip 
where one encounters the group of 5 Scuti stars defined 
here as Pop I stars in a postZams stage. The associated 
mass range approximately is 1.6-2 M Q . A mass of 1.80 
M© with a standard chemical composition for Population 
I (X=0.71, Y=0.27, Z= 0.02) is therefore chosen as typical 
for the present work. 

Fig.l shows the evolutionary track running from the 
latest pms stages to the early postZams ones for this 1.80 
M© stellar model. Peculiar stages are indicated and their 
characteristics listed in Table 1. Local maxima of the lu- 
minosity along the pms track correspond to the onset of 
rapid nuclear burning stages: ignition of the p-p reaction 
occurs close to models M6 and M7; the C 12 -N 14 transfor- 
mation (close to models M12, M13 and M14) gives rise to 
core convection. As the CNO cycle contribution becomes 
increasingly important and eventually dominant, the con- 
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Fig. 1. Evolutionary sequence of a 1.80 M Q model (solid 
line). Selected models listed in Table 1 are indicated by 
dots and labeled on the track. A 1.65 M© evolutionary 
sequence is also shown for comparison (dashed line). 



vective core slightly expands till the stage of M10. Once 
the M14 stage is reached, the star definitely settles into 
the CNO cycle burning phase. The zero age main sequence 
(Zams) model is defined as the model with the minimum 
luminosity just before evolution starts back toward the 
red. At this stage, M15 in Table 1, the gravitational en- 
ergy (e grav ) represents less than 1% of the total generated 
energy {ttot) and the central initial hydrogen content has 
only been slightly modified (see Table 1). 

This study does not include the earliest stages, hence 
it cannot provide exact values of ages. However, evolution- 
ary time scales for phases after the ignition of the central 
hydrogen burning phase are reliable. We have therefore 
arbitrarily chosen the stage Ml as the initial epoch. 

2.2. Structural differences between pms ans postZams 
models with same global properties 

As mentionned earlier, a source of ambiguity comes from 
the fact that pms and postZams stars have similar surface 
properties and therefore lie at the same position in a HR 
diagram. For instance, our 1.8M sequence in Fig.l shows 
that the three pms models M4, M8, M13 can be confused 
with the three postZams models M23, M19, M16 respec- 
tively, on the basis of their sole location in the HR diagram 
(Table 1). The question we address here is whether small 
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Table 1. Selected stellar models from the 1.80 M Q evo- 
lutionary sequence. The second column gives the effective 
temperature T of the model. Ages are in Myr, luminosi- 
ties L and total stellar radii, R, in solar units, convective 
radii r conv in unit of total stellar radii. X c is the central 
hydrogen content in percentage 
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log T 


log L 


R 


age 


x c 


fconv 


M 1 


3.808 


1.064 


2.751 


.0 


0.7143 


0.000 


M 2 


3.844 


1.133 


2.527 


.4 


0.7143 


0.000 


M 3 


3.863 


1.160 


2.389 


.6 


0.7142 


0.000 


M 4 


3.872 


1.170 


2.319 


.7 


0.7142 


0.068 


M 5 


3.880 


1.176 


2.251 


.8 


0.7142 


0.124 


M 6 


3.887 


1.177 


2.184 


.9 


0.7142 


0.166 


M 7 


3.895 


1.160 


2.056 


1.1 


0.7142 


0.221 


M 8 


3.897 


1.129 


1.970 


1.25 


0.7142 


0.243 


M 9 


3.890 


1.030 


1.816 


1.6 


0.7141 


0.258 


M10 


3.875 


0.924 


1.721 


2.4 


0.7139 


0.239 


Mil 


3.897 


1.003 


1.702 


3.8 


0.7136 


0.211 


M12 


3.912 


1.038 


1.657 


4.4 


0.7135 


0.186 


M13 


3.919 


1.050 


1.621 


4.8 


0.7134 


0.186 


M14 


3.923 


1.047 


1.586 


5.2 


0.7134 


0.201 


M15 


3.922 


1.006 


1.521 


11. 


0.7120 


0.205 


M16 


3.920 


1.044 


1.608 


300. 


0.6383 


0.204 


M17 


3.918 


1.058 


1.649 


400. 


0.6097 


0.202 


M18 


3.908 


1.102 


1.813 


700. 


0.5127 


0.197 


M19 


3.897 


1.131 


1.978 


900. 


0.4355 


0.191 


M20 


3.893 


1.139 


2.029 


950. 


0.4142 


0.190 


M21 


3.889 


1.146 


2.086 


1000. 


0.3920 


0.188 


M22 


3.876 


1.162 


2.248 


1120. 


0.3340 


0.182 


M23 


3.872 


1.167 


2.312 


1160. 


0.3132 


0.180 


M24 


3.866 


1.172 


2.383 


1200 


0.2913 


0.178 


M25 


3.851 


1.183 


2.594 


1300. 


0.2320 


0.170 
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Fig. 2. Behavior of the Brunt- Vaissala frequency N in 
normalized unit (see Eq.l) as a function of radius (nor- 
malized to the total radius) for models along the 1.80 M Q 
evolutionary sequence: N profiles for two pms models (Ml, 
M7), Zams model (M15) and two postZams models (M17 
and M25) are represented. 



differences in the structures of the pms and postZams 
models can induce significant changes in their oscillation 
spectra. The second issue is to determine whether these 
frequency differences can be used to discriminate between 
the two evolutionary stages. 

Accordingly, our procedure consists in comparing seis- 
mological properties of pms and postZams models having 
the same position in the HR diagram and the same mass 
(herafter 'associated models'). 

As expected, the outer layers of associated models are 
very similar: density and temperature profiles and there- 
fore the opacity profile are almost identical, we therefore 
refrain from showing them. 

On the other hand, the central regions significantly 
differ. A pms model is still slightly contracting and re- 
mains chemically homogeneous. In contrast, a postZams 
model has a core which traces its nuclear history, changes 
its chemical composition and density, and builds regions 
with (j, gradients. 

Our youngest models correspond to the phase of the 
beginning of the p — p chain, with a very small expanding 
convective core. During the main sequence phase, as the 
star ages, hydrogen is depleted in the central region and 



the fully mixed convective core receeds with time starting 
at the Zams stage (approximately M15). 

Signatures of these modifications are clearly seen on 
the profile of the Brunt Vaissala frequency, TV, as a func- 
tion of radius (Fig. 2). The sharp transition between zero 
values of N within the convective core and positive val- 
ues in the overlaying radiative layers delimits the edge of 
the convective core. As it can be seen in Fig. 2, during the 
pms evolution, the outward progression of the edge of the 
expanding convective core can be followed with the corre- 
lated outward progression of the inner edge of N. The in- 
ner N maximum also moves outward and decreases. When 
the model evolves beyond the Zams on the main sequence, 
the inner edge and maximum of the Brunt- Vaissala fre- 
quency receed inwards, while a narrower but higher peak 
associated with the fx gradient region forms. It can be ex- 
pected from Fig. 2 that modes with normalized frequencies 
(Eq.l below) in the range 4-10 are the most sensitive to 
the size of the convective core and consequently to the 
evolutionary status (pms or postZams) of the star. 
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3. Nonradial pulsations of pms and postZams 
models 

At every stage considered here, the time scale of evolution 
is much longer than the pulsation period. For instance our 
1.8M sequence indicates time intervals of about 0.1 — 1 
Myr between specific stages of evolution represented in 
Fig.l. Periods of oscillation of S Scuti stars typically range 
between lh- 3h, reaching roughly 20 h for unstable g- 
modes. Growth rates are of the order of 100-500 yr, still 
much smaller than evolutionary time scales. The oscilla- 
tion amplitudes are expected to grow or decay on a much 
faster (Kelvin- Hclmoltz) time scale than the (nuclear) evo- 
lutionary time-scale, so that the amplitudes have settled 
onto their saturating value (limit cycle) long before evo- 
lution can have any disturbing influence. These consider- 
ations authorize investigations of the pulsations of pms 
stars by means of usual linear nonadiabatic analyses of 



static model s at a given stage of evolution (Cox 1980, 
Unno ct al., 1989 ). We treat the nonadiabatic effects in 
the formulation of Unno et al. (1989), using Suran's os- 
cillation code. The oscillation eigenfrequencies (rad/s) are 
denoted as ui. We will also use frequencies in pHz, v, and 
normalized frequencies, <j defined as 

(1) 



y/GM/R 3 

Frequencies of the modes are labelled in a usual way: 
n,l, m for the radial order, degree and azimutal number. 
The radial order of a given mode is determined accord- 
ing to Christensen-Dalsgaard & Berthomieu (1992) with 
n = l for the fundamental radial mode. The growth 
rate is obtained directly from numerical computations 
but we rather consider the dimensionless r\ parameter 
(Stellingwerf, 1978; Dziembowski ( |l993| ):) 



IF 



77 = 



Jo \w(x)\dx 



(2) 



where x = r/R and W is the total mechanical work (Cox, 
1980; Saio & Cox, 1980) 
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w(x)dx 



7T 



M 



p u p ' 



" P 



)}dm (3) 



where SP, Sp are the Lagrangean variations of the pres- 
sure and density respectively. The work integral, W, is 
normalized with the total kinetic oscillation energy: 



E k = \u? f \Sr\ 2 p d 3 r 
2 Jo 



(4) 



where Sr is the Lagrangean radius variation. 

Interaction between pulsation and convection is ne- 
glected. This is certainly not justified for the cool stars 
but the effect is less important for stars close to the Zams. 
Neglecting this interaction equivalently affects pms and 
postZams models at the same effective temperature close 
to the Zams. Hence this cannot significantly affect the 
main conclusions of the present comparative study. 
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Fig. 3. Normalized frequencies of a 1.8 M model evolv- 
ing from stage 1 to stage 25 of Fig.l. Abcissa is R — R re f 
for postZams models and R re f — R for pms models (see 
text). For sake of clarity, I = 1 et I — 2 are plotted sep- 
arately: top frame I — 0, 1 and bottom frame £ = 0,2. 
Frequencies of radial modes with same radial order are 
connected with dotted lines; I = 1 and £ = 2 are con- 
nected with solid lines. Frequencies of associated pms and 
postZams models M4-M23; M8-M19; M13-M16 are pre- 
sented with open circles. 



3.1. Evolution of the power spectrum with stellar age 

Fig. 3 displays the evolution of the normalized frequencies, 
er, for a 1.8 M® model evolving from stage 1 (pms) to 
stage 25 (postZams) along the evolutionary track shown 
in Fig.l. The plot is restricted to I = — 2 modes. These 
modes are usually thought to be the easiest modes to be 
detected, as their visibility coefficients remain large af- 
ter integration over the whole stellar disk. The radial or- 
der interval is chosen such as to cover the range of vi- 
brationally unstable modes. Evolution is represented in 
abscissae by a quantity which is defined so as to change 
monotonously throughout the pms and main sequence 
stages: AR is the difference between the radius of the 
model, R, and the radius of a reference model, R e ff, which 
for obvious reasons is chosen to be the Zams model (M15 
model of Table 1). Time goes from left to right. For pms 
models, AR = R — R re f < and for postZams models 
AR = R ref - R > Q. 
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Normalized frequencies evolve with age quasi- 
symetrically with respect to the Zams stage. Frequencies 
of radial modes and nonradial pure p-modes with given n, t 
remain nearly constant when evolving from late pms to 
early postZams stages because the models change nearly 
homologously. 

On the other hand, g-modes, as they are sensitive to 
the structure of the inner convective layers, do not behave 
homologously. One can see in Fig. 3 changes of the frequen- 
cies of the g-modes which arise from the 'breathing' of the 
core, (which moves back and forth) when a new nuclear 
reaction starts, at models M4 and M12: frequencies of g- 
modes increase with age for postZams models because the 
inner maximum of the Brunt- Vaissala frequency increases 
with age and its sharp inner edge moves outward. In a 
symmetric way with respect to the Zams model, the fre- 
quencies of (/-modes decrease when the pms model evolves 
toward the Zams stage because the position of the inner 
maximum of N moves inward with time. 

In addition, the well known phenomenon of avoided 
crossing (Unno et al, 1989) occurs during the latest 
postZams stages. The development of a \i gradient at the 
edge of the convective core of the postZams model causes 
the development of a sharp peak of the Brunt- Vaissala fre- 
quency Then nonradial low frequency modes undergo an 
avoided crossing which shifts the frequency with respect to 
its value in absence of avoided crossing. Only those modes 
which are sensitive to the edge of the convective core have 
a different behavior before and after the Zams stage. These 
interacting modes exchange their physical nature. Their 
frequencies penetrate into the p mode frequency domain, 
breaking the regularity seen at high frequencies. The pro- 
gression of this phenomenon into the p-mode domain can 
be followed from mode with n ~ 1 (AR ~ 0) up to n = 8 
(AR ~ 1) for I = 2 modes. This leads to different and 
recognizable patterns in a frequency spectrum, which can 
be considered as the signature of different structures of 
the core and can therefore tell whether the star is a pms 
or a main sequence one. 

Finally, frequency calculations show the existence of 
strange modes i.e. surface modes trapped in very superfi- 
cial outer layers (Unno et al. 1989; Buchler et al. 1997; Saio 
et al. 1998). Whether these modes are excited or not is 
probably strongly dependent on the convection-pulsation 
interaction not included here. As these surface modes are 
not essential for our purpose, they were disregarded in 
Fig.3. 



3.2. Vibrational instability 

The linear nonadiabatic calculations for modes with t = 
0,1,2 for models of Table 1 show that the pms models 
lying in the classical instability strip are vibrationally un- 
stable. The unstable modes are the same than those for 
classical 5 Scuti stars, namely low radial order g- and p- 
modes. 
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Fig. 4. Growth rate r\ versus normalized frequency a for 
low radial order I = 0,1,2 modes for the associated pms 
(M8, open circles) and main sequence (M19, full squares) 
models. 



Fig. 4 shows the behavior of the parameter r\ (Eq.2) 
with increasing normalized frequencies a for I = 0,1,2 
modes for the M8 pms and the M19 main sequence models. 
The instability strength as measured with 77 is independent 
of the degree I. The 77 curves for pms and main sequence 
associated models are almost identical. For these models, 
modes with radial order in the range n = 1 — 7 are found 
to be unstable. 

The work integrands are found to be quite similar for 
associated pms and postZams models. Computation of the 
work integrals shows that, for both pms and postZams 
models, the modes are driven by opacity mechanisms 
which efficiently operate in the H and He ionisation zones. 
The e mechanism contributes only very weakly in the deep 
interior even in the youngest models. 

4. Frequency distributions of pms versus 
postZams models 

Fig. 5 compares frequencies of modes with the same degree 
i for the associated models M4 and M23, represented re- 
spectively by solid vertical lines and dashed vertical lines. 
Frequencies of I — modes with same radial order are 
nearly identical as shown by the solid lines which nearly 
coincide with the dashed lines. This is because both mod- 
els have similar mean density and outer layers. This is also 
true for the highest frequencies of I = 1 and I = 2 modes 
which are associated with pure p modes. Differences in 
the spectral patterns on the other hand exist in the low 
frequency range. In particular a few extra frequencies, as- 
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Fig. 5. Frequency distribution for £ = 0,1,2 modes for 
pms M4 (solid lines) and main sequence M23 (dashed 
lines) models. Ordinate is arbitrary. Abscissa is the nor- 
malized frequency a. 

sociated with £ = 1,2 modes in avoided crossing, exist 
for the postZams model. As a consequence, even without 
a mode identification, a comparison of the spectra for a 
pms and a postZams models at the same location in the 
HR diagram tells us whether the star is in a pms or a 
postZams stage and, in the latter case, yields an indica- 
tion of its age. 

Frequency differences for a given mode between mod- 
els M4 and M23 reach values as large as 50 /iHz in real 
units. This is emphasized in Fig. 6. Global parameters are 
identical, therefore radial modes with same radial order 
have about the same frequencies. The external layers are 
similar, therefore the excitation range is about the same. 
In this frequency range, £ = 2 modes are the most affected 
by the structure of the deep layers and are the main agents 
for differences in power spectra. 

5. Toward a seismological diagnostic: frequency 
differences 

Pms and postZams models at the same location in the 
HR diagram have similar mean densities, hence one ex- 
pects the large separation v n ^ — v n -i,t to be the same for 
'associated' models. This is confirmed in Fig. 7 where the 
large separations for t = modes (for low radial order 
modes) are displayed for associated models and for the 
Zams model. The large separation is indeed the same re- 
spectively for M4 and M23; for M8 and M19; for M13 
and M16. The large separation scales as (GM/R 3 ) 1 / 2 . 
Given the mass, the large separation increases with age 
for pms models (decreasing radii) and decreases with age 
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Fig. 6. Theoretical power spectra in the low frequency 
range where modes are excited, top: postZams M23 
model, bottom: pms M4 model. Modes with same (£, n) 
subscripts in both models are connected with dashed lines 
(£ = 0,£ = 1 modes) and dotted lines (£ = 2) modes. 
Amplitudes are arbitrary. For clarity the highest ampli- 
tudes correspond to I = modes (thicker lines), medium 
ones to I = 1 modes and smallest ones to I = 2 modes. 
Marginally stable modes (modes with —0.005 < r\ < 0) 
are pictured with dash-dotted vertical lines. 

for postZams models (increasing radii); it is largest for the 
Zams model which has the smallest radius. 

The development of a (i gradient at the edge of the 
convective core can be seen on the frequency large sepa- 
ration for nonradial modes in Fig. 8. The large separation 
for 1 = 1 modes for associated pms and postZams models 
is expressed in normalized units and plotted against the 
normalized frequency (J n ,e=2- The existence of a mode in 
avoided crossing generates a large departure of the fre- 
quency difference a n ^ — a n -i,e from its mean value at the 
frequency of this mode. The frequency at which a given 
mixed mode occurs is correlated with the maximum value 
of the inner peak of the Brunt- Vaissala, frequency arising 
from a larger [i gradient. Such a mixed mode is detected 
in postZams model M20, M22, M23 and M24 in Fig.8. 
While the main sequence model ages, the progression of 
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Fig. 7. Large separations v n e — Vn-i,i ip-Hz) for low ra- 
dial order I = modes for associated models from bot- 
tom to top M4-M23; M8-M19; M13-M16 (pms:solid lines; 
postZams dashed lines) and the Zams model Ml 5 (dotted 
line) : large separations are the same for pms and postZams 
models with the same mean density. 



the mixed mode toward higher frequencies is clearly seen. 
The frequency of the extra mode with the highest fre- 
quency gives an indication of the age of the postZams star. 
On the other hand, pms models show a quasi constant 
large separation (in normalized units) with frequency. 

Differences in the inner layers of pms and postZams 
models can also be detected by investigating the behav- 
ior of the small separations v n j — v n ',e-2 as a function of 
the frequency v n j. In the asymptotic regime i.e. for high 
radial order modes (solar like oscillations) the small sepa- 
ration is defined as v n ^ — i/ n _i^_ 2 . In the low radial order 
range, avoided crossings perturb the regular sequence of 
alternations of I = 2 and 1 = modes. We therefore find 
more convenient to define a small separation as the dif- 
ference between a I = 2 mode (radial order n') and its 
closest t = neighbor (radial order n which in some cases 
at low frequency differs from n' + 1). Such small separa- 
tions are plotted in Fig. 9. When the associated models are 
very close to the Zams as for instance models M13-M16, 
no difference exists between the small separations of the 
models. For these models, as no \i gradient has developed 
yet, there is no avoided crossing in the excited frequency 
range. 

On the other hand, for associated models located fur- 
ther away from the Zams, for instance M4-M23, the struc- 
tures of the inner layers arc significantly different; this 
leads to a different behavior of the small separations of the 
two models at low frequency. An avoided crossing is seen 



Fig. 8. Large separations a n .g — cr n -i,i f° r (■ — 2 ver- 
sus (T n ,e=2 for associated pms (dotted lines) and postZams 
(solid lines) models from top to bottom: M13-M16; M7- 
M20; M5-M22; M4-M23; M3-M24. 

to occur for modes with v ~ 350^iHz [a ~ 8) for model 
M23 which is not present for its associated pms model M4 
(Fig. 9). At lower frequencies (below 250^Hz correspond- 
ing to cr ~ 6), a second feature is seen for both models 
M4 and M23: this sharp variation of the small separation 
is due to modes with frequencies close to the value of the 
secondary inner maximum of the Brunt- Vaissala frequency 
(the only inner maximum for pms models) near the edge 
of the convective core. 

6. Discussion 

Pre-main sequence models of 1.7-2 M Q crossing the classi- 
cal instability strip on their way to the Zams, are found to 
be vibrationally unstable, destabilized by the same mech- 
anisms than the S Scuti stars, meant as main sequence or 
postZams stars. 

For a given mass and a given location in the HR di- 
agram, the frequencies of radial modes with same radial 
order remain nearly identical. As the outer layers in both 
stages are almost the same, oscillations in pms stars should 
have the same level of amplitudes as in 8 Scuti stars; this 
means that some modes are expected to be detectable 
from ground. 

Low radial order modes with dual (gravity and pres- 
sure) nature penetrate in the upper layers of the star and 
there is a good chance to observe them. We have shown 
that they can act as discriminators, as they are very sen- 
sitive to the structure of the inner layers, which differ be- 
tween pms and postZams stars. This is particularly inter- 
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Fig. 9. Small separations v n ,i=o — Vn',t=2 versus v n ^ = 2 
(in [iHz) for associated pms (solid lines, open circles) 
and postZams (dashed lines, open stars) models: M4-M23 
(bottom) and M13-M16 (top). The frequency difference is 
calculated between a I = 2 mode and the I = mode with 
the closest frequency. 

esting for this region of the HR diagram where classical 
indicators of youth as accretion disks or remaining nebu- 
losities have disappeared. 

Frequency shifts and differences in the frequency dis- 
tributions of pms and main sequence oscillators can be 
detected if a sufficient number of modes are observed at 
low frequencies in the vinicity of the fundamental radial 
mode. 

Important progress has been made to detect larger sets 
of frequencies from the ground using networks of tele- 
scopes covering all longitudes. Nevertheless ideal condi- 
tions will only be achieved in space. The COROT exper- 
iment (Baglin et al. 1998), to be launched in 2004, will 



continuously observe several stars during 20 to 150 days; 
this will provide an increase of a thousand for the signal 
to noise ratio and a frequency resolution from 0.5 fj,Hz 
down to 0.1 [iHz . The large number of modes which is 
expected to be detected will help to provide, at least, a 
partial mode identification. Frequency histograms for in- 
stance can provide the mean density (Breger et al, 1999; 
Baglin et al, 2000; Barban et al. 2001) and the structure 
of the frequency spectrum will tell whether the star is a 
pms or a postZams star. 

Pms stars are expected to be fast rotators and it will 
be necessary to include effects of rotation on the frequen- 
cies in the same way as proposed by Soufi et al. ( 1998 ) 
for fast rotating 5 Scuti stars. This should enable to put 



constraints on the internal rotation of these stars (Goupil 
et al., 1996| ). It will then be possible to trace the history- 
and to understand better the process -of transfer of an- 
gular momentum inside stars with ~ 2Mq in their early 
stages of evolution. 
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